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RMS  LG  MEASUREMENTS  OF  CHINESE  UNDERGROUND 
EXPLOSIONS  AT  LOP  NOR 


Won- Young  Kim  and  Paul  G.  Richards* 
Lamont-Doherty  Earth  Observatory,  Palisades,  NY  10964 
(*  also.  Dept,  of  Geological  Sciences,  Columbia  University) 


ABSTRACT 

RMS  Lg  measurements  using  digital  seismograms  recorded  at  the  Borovoye  Geophysical 
Observatory,  Kazakhstan  (station  BRVK)  from  Chinese  underground  nuclear  explosions  in  the 
Lop  Nor  test  site  (A  «  1880  km)  show  fairly  good  correlation  with  ISC  m^(P).  Regression  of 
five  RMS  Lg  measurements  using  BRVK  data  yields  a  standard  deviation  of  only  0.045 
magnitude  units. 


INTRODUCTION 

In  this  study,  we  used  vertical  component  digital  seismograms  from  five  Lop  Nor 
underground  explosions  (magnitude  4.7  -  6.2)  recorded  on  short-period  seismographs  at  BRVK. 
Epicentral  distances  ranges  from  1879  to  1894  km  with  a  mean  distance  of  1884  km  (Table  1). 
Four  explosions  records  on  SKM-3  high-gain  channel  of  STsR-SS  seismograph  and  an  explosion 
record  on  KSVM  low-gain  channel  of  STsR-TSG  seismograph  are  used  for  RMS  Lg 
measurements.  Characteristics  of  the  STsR-SS  and  -TSG  seismographs  are  given  in  Table  2. 

The  SS-SKM-3  channels  use  Kimos  short-period  seismometers  and  record  three-component 
high-gain  and  vertical  component  low-gain  channels.  The  high-gain  channels  have  been  operating 
with  a  nominal  gain  of  about  2000  count/p.  and  have  nearly  flat  response  to  ground  displacement  in 
the  frequency  band  0.8  -  3.3  Hz  (3  dB  level).  This  high-gain  vertical  channel  provides  good 
regional  signals  with  frequencies  up  to  about  3-4  Hz  for  most  of  the  underground  explosions  from 
the  Lop  Nor  test  site  (Fig.  1).  The  KSVM  channel  of  the  TSG  seismograph  has  similar 
characteristics  as  SS-SKM-3  channels,  but  has  only  low-  and  high-gain  vertical  components. 


METHOD 

To  measure  the  RMS  Lg,  we  followed  a  procedure  similar  to  Hansen  et  al.  (1990),  except 
that  we  used  a  gaussian  window  instead  of  the  usual  box-car  window.  Noise  in  the  signal  is 


corrected  as  suggested  in  Ring  dal  &  Hokland  (1987)  by  taking  RMS  values  of  the  trace  (50  sec 
window)  preceding  the  first  arrival  P-wave.  Lg  signal-to-noise  ratios  were  about  2.5  to  92. 

We  used  a  gaussian  with  width  cStef  =  45  s  at  a  reference  distance,  Aref  =  1000  km.  This 
gaussian  window  (icr)  covers  the  Lg  phase  in  group  velocity  range  3.66  to  3.0  km/s  when 
centered  at  the  group  velocity  around  3.3  km/s.  The  gaussian  is  truncated  at  2.58a  (99  %  of  unit 
area).  The  width  of  the  gaussian  can  be  conveniently  scaled  as  a  function  of  distance,  a  =  Gref  * 
A/Aref,  allowing  the  RMS  Lg  measurements  at  many  stations  at  different  epicentral  distance  ranges 
to  be  combined  to  obtain  a  network  m^(Lg)  based  on  the  RMS  Lg  measurements.  The 
successive  gaussian  windows  were  shifted  by  0.2  a  relative  to  the  previous  window  to  sample  the 
Lg  waves  smoothly.  The  Lg  signals  are  not  corrected  for  the  geometrical  spreading  and  anelastic 
attenuation  along  the  path  because  all  explosions  are  at  nearly  the  same  location. 


RESULTS 

Fig.  2  shows  the  logio  of  RMS  Lg  measurements  in  nanometer  of  ground  displacement 
plotted  against  ISC  nti,(P).  Regression  of  ISC  m^(P)  against  logio  (RMS  Lg)  for  five 
explosions  in  the  magnitude  range  4.7  -  6.2  yields  a  slope  of  0.98  and  a  standard  deviation  of  only 
0.045  magnitude  units.  It  appears  that  BRVK  data  is  useful  to  measure  source  strength  of  Lop 
Nor  explosions  for  explosions  with  magnitude  down  to  about  m^(P)  =  4.5.  Examination  of  all 
available  GSN  (Global  Seismograph  Network)  digital  seismograms  for  13  known  underground 
explosions  during  1976-1992  at  Lop  Nor  indicates  that  BRVK  provides  unique  regional  data  for 
RMS  Lg  studies. 


Table  1.  Lop  Nor  Underground  Explosions(*) 


Origin  time 

Latitude  Longitude 

ISC  Magnitude 

logio(RMS  Lg) 

S/N 

Year  Mon  Day  ( h:  m:  sec) 

(°N) 

(°E) 

("lb) 

# 

(Ms) 

# 

(nanometer) 

1983  OCT  06 

09:59:58.0 

41.53 

88.72 

5.5 

73 

4.2 

2 

1.3907 

13.7 

1984  OCT  03 

05:59:57.9 

41.54 

88.67 

5.4 

81 

- 

1.1956 

10.8 

1984  DEC  19 

06:00:02.8 

41.62 

88.22 

4.7 

11 

4.2 

1 

0.6236 

2.5 

1987  JUN  05 

04:59:58.5 

41.55 

88.72 

6.2 

145 

4.7 

5 

2.1385 

92.3 

1990  AUG  16 

04:59:57.7 

41.52 

88.75 

6.2 

138 

4.4 

17 

2.1008 

67.4 

(*)  Origin  time,  location  and  magnitudes  are  from  ISC  bulletin;  S/N=signal-to-noise  ratio. 


Vertical-component  Records  at  BRVK  from  Lop  Nor  Underground  Explosions 


Fig.  1.  Vertical  component  records  (SS-SKM-3)  recorded  at  BRVK  from  Lop  Nor  explosions. 


The  RMS  measurements  of  P  waves  are  also  obtained  following  a  similar  procedure  used 
for  Lg  signals  to  examine  the  stability  of  the  method  for  the  P  signals.  Regression  of  ISC 
mfjiP)  against  logio  (RMS  P)  values  at  BRVK  for  six  explosions  in  the  magnitude  range  4.7  - 
6.2  yields  a  slope  of  1.03  and  a  standard  deviation  of  0.065  m.  u.  Though  this  result  is 
preliminary  in  nature,  it  indicates  that  RMS  P  wave  signals  may  provide  an  alternative  tool  for 
measuring  the  strength  of  the  seismic  source  when  the  records  have  only  P-wave  portion  of 
triggered  data. 


Table  2.  Characteristics  of  STsR  Seismographs  at  Borovoye  Seismic  Station^*). 


Seismograph  Seismometer  Data 

np  (a) 

As 

c  (c) 

f  (d) 

im 

dt^*^) 

Channel® 

channel 

(s) 

(unit/|X) 

(Hz) 

(msec) 

number 

STsR-SS 

SKM-3 

High-gain 

2.0 

0.5 

2000 

0.8-3.3 

24 

7,8,9 

Low-gain  (Z)t 

200 

32 

1 

20 

96 

6 

SKD 

High-gain 

25.0 

0.5 

5.0 

0.05-0.4 

192 

2,3,4 

Low-gain 

0.5 

192 

1,5,10 

STsR-TSG 

KS 

1.5 

0.7 

4500 

1.0-3.3 

26 

7,8,9 

DS 

20.0 

0.5 

50 

0.05-0.2 

312 

19,20,21 

KSM 

High-gain 

1.5 

0.35 

100000 

0.7-5.0 

26 

10,11,12 

Low-gain 

1000 

26 

3,4,5 

DSM 

High-gain 

28.0 

0.5 

1000 

0.04-0.1 

312 

22,23,24 

Low-gain 

10 

312 

15,16,17 

KSVM 

High-gain  (Z) 

1.5 

0.35 

4600 

0.8-2.5 

26 

2 

Low-gain  (Z) 

50 

26 

1 

(*)  STsR  system  from  Feb  1973  to  present,  (a)  Ts=Seismometer  natural  period  in  second,  (b) 
Ds=Seismometer  damping  constant,  (c)  Sjii=Nominal  gain  in  count/|i,  (d)  fni=passband  (90  %  of 
peak),  (e)  dt=Sampling  interval  in  millisecond,  (e)  Channel  identifier,  t  Only  vertical  component . 
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Fig,  2  Comparison  of  logjo  of  RMS  Lg  (in  nm)  at  BRVK  with  ISC  m^CP).  Solid  line  is  a  fitted 
slope  of  0.98  and  an  orthogonal  rms  misfit  of  0.045  m.  u.  Dotted  lines  correspond  to  ±2  S.D. 


CONCLUSIONS 

The  regional  seismic  data  now  becoming  available  for  Central  Asia,  digitally  recorded  at 
Borovoye,  Kazakhstan,  appear  to  be  useful  to  measure  source  strength  of  Lop  Nor  exf  losions  for 
explosions  with  magnitude  down  to  about  »  4.5. 

Instrument  calibration  (gains)  of  various  seismographs  at  Borovoye  station  appear  to  be 
stable  over  many  years.  Stability  of  Lg  amplitudes  for  explosions  at  the  Lop  Nor  test  site  over 
several  years  suggests  reasonable  instrument  calibration.  The  broadband  nature  of  BRVK  digital 
records  provides  an  opportunity  to  examine  seismic  velocities  and  discontinuities  in  the  upper 
mantle  beneath  Central  Asia.  In  these  conclusions,  the  common  theme  is  new  opportunities  to 
work  with  high-quality  regional  signals  for  Central  Asia  recorded  at  BRVK.  We  recommend  a 
concerted  effort  by  academic  institutions  to  ensure  timely  salvage  of  these  valuable  seismic  data. 
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INSTRUMENT  RESPONSES  OF  DIGITAL  SEISMOGRAPHS  AT 
BOROVOYE,  KAZAKHSTAN  BY  INVERSION  OF  TRANSIENT 

CALIBRATION  PULSES 


Won- Young  Kim^  and  Goran  Ekstrom^ 
Lamont-Doherty  Earth  Observatory  of  Columbia  University, 
Department  of  Earth  and  Planetary  Sciences,  Harvard  University 


EXTENDED  ABSTRACT 

A  method  is  developed  to  determine  the  response  of  digital  seismographs  from  their  transient 
calibration  pulses  when  there  is  no  recourse  to  direct  measurements.  Based  on  linear  system 
theory  for  small  input  signal,  the  digital  seismograph  is  represented  by  a  set  of  first-  and  higher- 
order  linear  filters  characterized  by  their  cutoff  frequency  and  damping  constant.  The  transient 
calibration  pulse  is  parameterized  by  a  set  of  instrument  constants,  and  the  problem  is  linearized  for 
small  perturbation  with  respect  to  their  nominal  values.  The  calibration  pulse  shape  is  matched  in 
the  time  domain  using  an  iterative  linearized  inverse  technique. 

The  method  is  used  to  derive  complete  instrument  responses  for  digital  seismographs 
operating  at  the  Borovoye  (BRVK)  in  Kazakhstan,  for  which  previously  only  the  amplitude 
responses  have  been  determined.  To  test  the  method,  we  apply  it  to  digital  calibration  pulses  from 
a  modem  digital  seismograph  system  at  Kislovodsk  (KIV)  in  northern  Caucasus,  Russia,  and 
obtained  good  agreement  between  known  and  derived  instrument  constants.  The  results  of  the 
calibration  pulse  shape  inversion  for  these  seismographs  indicates  that  the  method  is  efficient  and 
that  the  results  are  reliable  even  when  microseismic  noise  is  present  in  the  recorded  transient 
calibration  pulse.  The  derived  parameters  make  possible  improved  quantitative  waveform  analysis 
of  digital  seismograms  recorded  at  BRVK. 

In  this  extended  abstract,  we  present  only  summary  tables  (Table  1  &  Table  2)  and  a  figure 
for  the  STsR-TSG  seismic  system  which  operated  since  1973  at  BRVK.  The  full  text  with  figures 
and  tables  will  appear  in  February  1996  issue  of  the  Bulletin  of  the  Seismological  Society  of 
America. 


Table  1.  Characteristics  of  STsR-TSG  Seismic  System  at  Borovoye  Seismic  Station^*). 


Seismograph  Data 

channel 

yCa) 

is 

(S) 

5s® 

Gain^^^ 

(count/p) 

f  (d) 

In 

(Hz) 

dt<'> 

(msec) 

Channel^^ 

number 

Valid® 

date 

KS 

1.5 

0.7 

2250 

1.5 

26 

7,8,9 

07/23/76  - 

1.5 

0.7 

4500 

1.5 

26 

7,8,9 

03/24/82  - 

DS 

20.0 

0.5 

50 

0.1 

312 

19,20,21 

07/12/74  - 

KSM 

hg(*'> 

1.5 

0.35 

100000 

1.0 

26 

10,11,12 

11/01/81  - 

LG® 

1000 

26 

3,4,5 

DSM 

HG 

28.0 

0.5 

1000 

0.07 

312 

22,23,24 

09/08/82  - 

LG 

10 

312 

15,16,17 

KSVM 

HG  (Z)t 

1.5 

0.35 

4600 

1.0 

26 

2 

12/12/83  - 

LG  (Z)t 

50 

26 

1 

(*)  STsR-TSG  system  from  Feb  1973  to  present,  (a)  Ts  =  Seismometer  natural  period  in  second, 
(b)  =  Seismometer  damping  constant,  (c)  Gain  =  Nominal  gain  in  counts/p,,  (d)  f^  = 

normalization  frequency  where  the  nominal  gain  is  measured,  (e)  dt  =  Sampling  interval  in 
millisecond,  (f)  Channel  identifier  on  the  original  station  tape,  (g)  dates  when  the  responses  given 
are  valid,  (h)  High-gain  channel,  (i)  Low-gain  channel,  f  Only  vertical-component 


Table  2.  Instrument  constants  of  the  STsR-TSG  seismographs  from  inversion*. 
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Table  2.  Continued. 
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and  the  nominal  gain  since  03/24/82  -  present  is  4,500  counts/^t; 

(b)  KSM  low-gain  channels  have  identical  responses,  but  with  lower  gains  of  Z  (3)  =  982,  NS  (4)  =  1002,  EW  (5)  =  995; 

(c)  KSVM  low-gain  channel  has  identical  response,  but  with  a  lower  gain  of,  Z  (2)  =  46.15; 

(<l)  DSM  low-gain  channels  have  identical  responses,  but  with  lower  gains  of,  Z  (15)  =  10.0,  NS  (16)  =  9.9,  EW  (16)  =  10.0. 


STsR-TSG  Seismic  System 


Frequency  (Hz) 

Fig.  1.  A  summary  figure  showing  spectral  amplitude  responses  of  the  all  vertical-component 
seismographs  of  the  TSG  system.  In  each  seismograph,  frequency-amplitude  responses  given  in 
the  log  book  at  BRVK  {closed  circles)  and  the  amplitude  responses  obtained  after  the  calibration 
pulse  waveform  inversion  {solid  lines)  are  plotted  together  for  comparison. 
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ABSTRACT 

Several  Russian  sources  have  stated  that  343  underground  nuclear 
explosions  were  conducted  during  1961-1989  at  the  Semipalatinsk  Test 
Site.  However,  only  282  of  them  appear  to  have  been  described,  in  the 
openly  available  technical  literature,  with  well-determined  coordinates; 
and  only  272  have  both  good  locations  and  magnitudes. 

We  have  used  regional  data  from  52  stations  to  study  65  seismic 
sources  initially  thought  to  be  in  or  near  the  Semipalatinsk  region, 
additional  to  the  272  underground  nuclear  explosions  with  known 
locations  and  magnitudes.  Of  these  65  events,  we  believe  8  are  not 
explosions  on  the  test  site,  namely;  two  earthquakes  close  to  the  test  site 
(one  of  them,  on  March  20,  1976,  was  previously  well-known);  three 
earthquakes  or  chemical  explosions  100-300  km  from  the  test  site;  and 
three  events  at  greater  distances  from  Semipalatinsk.  Of  the  remaining  57 
events:  10  were  known  to  be  underground  nuclear  explosions  with  known 
locations  and  we  have  supplied  magnitudes  where  none  were  previously 
available;  one  was  a  chemical  explosion  at  Degelen  (June  5,  1961,  a  few 
months  prior  to  the  first  Soviet  underground  nuclear  explosion);  we 
believe  21  were  underground  nuclear  explosions  (20  at  Degelen,  one  at 
Murzhik);  13  were  chemical  explosions  at  Balapan;  8  were  chemical 
explosions  elsewhere  on  the  test  site;  three  were  either  nuclear  or  chemical 
explosions;  and  one  was  either  a  chemical  explosion  or  a  cavity  collapse. 

Our  seismological  data  is  principally  of  two  types;  (1)  the  bulletins  of 
stations  in  Central  Asia,  Kazakhstan  and  the  Altai;  and  (2)  the  multi¬ 
channel  narrowband  analog  records  (ChISS)  of  the  station  at  Talgar, 
operated  throughout  the  time  period  1961-89,  together  with  other  ChISS 
stations  at  Garm,  Zerenda,  and  Novosibirsk,  operated  for  part  of  this  time. 
The  largest  magnitude  of  our  44  possible  underground  nuclear  explosions 
is  around  5  (February  4,  1965,  obscured  at  many  teleseismic  stations  by  a 
large  Aleutian  earthquake).  Others  lie  in  the  magnitude  range  3.5-4. 5,  and 
clearly  most  have  subkiloton  yields. 
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Our  data  set  of  small  events  is  important  for  purposes  of  evaluating 
the  detection  capability  of  teleseismic  arrays,  and  the  detection  and 
identification  capability  of  regional  stations. 


INTRODUCTION 

This  paper  reports  our  practical  experience  in  documenting  the 
occurrence  of  small  underground  nuclear  explosions  and  large  chemical 
explosions  at  the  main  nuclear  weapons  test  site  of  the  U.S.S.R.,  near 
Semipalatinsk,  Kazakhstan. 

We  have  worked  with  five  sources  of  information.  First,  are  the 
various  lists  of  explosions  published  by  western  scientists  based  upon 
teleseismic  observations.  Second,  are  lists  of  Soviet  explosions  and 
summary  statements  on  total  numbers  of  such  explosions,  published  by 
Russian  sources.  We  find  that  teleseismically  observed  events  appear  to 
be  only  about  90%  of  the  total.  Some  of  the  teleseismically  observed 
events  are  poorly  located,  or  have  not  been  assigned  a  magnitude,  or 
appear  on  explosion  lists  that  incorrectly  include  earthquakes.  It  appears 
that  less  than  80%  of  the  underground  nuclear  explosions  conducted  at  the 
Semipalatinsk  Test  Site  have  good  locations  and  a  magnitude,  determined 
and  published  on  the  basis  of  teleseismic  data.  Third,  we  have  used  our 
experience  in  the  analysis  of  data  recorded  at  the  Talgar  station  on  a  day- 
to-day  basis,  and  at  stations  on  the  Kokchetav  massif  (Zerenda,  Borovoye). 
Fourth,  we  used  regional  seismic  data  from  49  other  stations  in 
Kazakhstan,  Kyrgystan,  Russia,  Tadjikistan  and  Uzbekistan,  most  of  them 
operated  at  distances  of  less  than  1000  km  from  the  Semipalatinsk  Test 
Site  in  East  Kazakhstan.  Fifth,  we  have  used  personal  communications 
from  western  seismologists  (Frode  Ringdal,  Peter  Marshall,  Jack  Murphy, 
Bob  North)  giving  us  information  based  for  the  most  part  on  teleseismic 
data  from  particular  events. 

Our  goals  were  to  reduce  the  discrepancies  between  explosion  lists 
generated  from  teleseismic  data  and  numbers  of  events  reported  by 
Russian  sources,  and  thus  to  develop  a  more  complete  list  of  Soviet  nuclear 
explosions;  and  to  document  some  strengths  and  weaknesses  of  teleseismic 
and  regional  data  for  purposes  of  explosion  detection  and  identification. 
More  fundamentally,  we  wanted  to  build  a  database  of  small  seismic 
events,  so  that  efforts  can  be  focussed  on  solving  the  hardest  technical 
problems  of  monitoring  a  Comprehensive  Test  Ban  Treaty,  namely  how  to 
identify  underground  nuclear  explosions  that  have  very  small  seismic 
signals  -  and  how  to  avoid  false  alarms  over  small  events.  This  work  is 
complicated  in  the  case  of  seismic  events  near  Semipalatinsk,  by  the  fact 
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that  75  chemical  explosions  have  been  carried  out  on  the  test  site, 
according  to  Saybekov  (1993),  including  44  chemical  explosions  larger  than 
ten  tons  that  have  been  conducted  during  the  period  1970-1988. 

Our  database  so  far,  as  reported  in  this  paper,  consists  of  numerous 
station  reports  and  analog  seismograms  that  we  have  analysed  ourselves. 
We  often  relied  upon  methods  of  discrimination  that  cannot  in  general  be 
used  in  CTBT  monitoring  (for  example,  the  interpretation  of  the  time  of  day 
at  which  the  event  occurred).  However,  we  believe  we  have  identified  a 
suitable  set  of  small  events  (both  explosions  and  earthquakes),  on  and  near 
the  Semipalatinsk  Test  Site,  that  may  be  an  important  focus  for  future 
work.  That  is,  we  anticipate  further  efforts  to  acquire  data  for  these 
events,  preferably  digital  data,  recorded  at  teleseismic  and  regional 
distances,  in  order  to  test  objectively  the  capability  of  various  detection 
and  discrimination  algorithms.  In  such  future  work  we  could  digitize  some 
of  the  analog  data  reported  here,  or  we  could  hope  to  use  digitally 
recorded  data  from  Borovoye  since  1965  (at  a  distance  about  700  km), 
some  digital  data  acquired  by  the  Natural  Resources  Defense  Council  in 
1987-1988  (at  distances  about  200  km  from  the  test  site),  and,  since  1988, 
IRIS/IDA  digital  data  (at  distances  more  than  1000  km). 


SOURCES  OF  INFORMATION  ABOUT  SMALL  EVENTS 

In  this  section,  we  first  summarize  the  information  contained  in  a 
series  of  papers  published  in  the  West,  concerning  numbers  of 
underground  nuclear  explosions  (UNEs)  that  have  occurred  at  the 
Semipalatinsk  Test  Site  in  Kazakhstan,  and  for  which  accurate  locations  and 
magnitudes  are  both  available.  We  then  briefly  summarize  additional 
information  made  available  in  recent  years  by  various  Russian  and  Kazakh 
sources,  which  indicates  that  a  significant  number  of  UNEs  at  Semipalatinsk 
have  not  been  previously  recognized  in  western  publications.  The  section 
concludes  with  a  list  of  65  events  which  are  the  main  focus  of  attention  in 
this  paper,  being  candidates  for  the  additional  UNEs,  and/or  UNEs  with 
previously  unknown  locations  or  magnitudes,  or  perhaps  earthquakes  or 
chemical  explosions  in  the  region. 

Table  1  gives  the  information  available  on  282  known  UNEs  at  the 
Semipalatinsk  Test  Site.  Our  primary  sources  of  information  on  locations 
comes  (1)  from  Bocharov  et  al  (1989),  who  published  origin  time,  depth, 
latitude  and  longitude  on  96  UNEs  prior  to  1973;  and  (2)  from  Lilwall  and 
Farthing  (1990)  who  used  the  Bocharov  et  al  locations  as  master  events  for 
a  joint  epicentral  determination  of  other  UNEs  at  Semipalatinsk,  using  P- 
wave  arrival  times  as  reported  by  the  ISC.  The  Lilwall  and  Farthing 
seismically  determined  locations  were  given  to  the  nearest  thousandth  of  a 
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degree  (five  significant  figures).  The  "area"  entry  in  Table  1  is  the  area  of 
the  error  ellipse  given  by  Lilwall  and  Farthing  (blank  entries  signify 
events  given  a  location  to  the  nearest  tenth  of  a  second  of  arc  by  Bocharov 
et  al:  such  locations  are  listed  in  Table  1  with  seven  significant  figures). 

An  event  of  1984  September  15,  presumed  to  be  a  UNE  located  at 
49.992  °N,  78.881  °E  by  Lilwall  and  Farthing,  was  apparently  a  chemical 
explosion  at  Balapan.  The  event  of  1968  November  12  was  a  triple 
explosion  (Bocharov  et  al)  so  we  have  counted  it  three  times. 

Table  1  also  indicates  the  subsite  (B  =  Balapan,  D  =  Degelen,  M  = 
Murzhik);  the  ISC  F-wave  magnitude,  where  available;  the  number  of 
stations  upon  which  the  ISC  magnitude  is  based;  the  number  of  stations 
reporting  dilatations  and  the  number  reporting  polarity;  and  the  number 
of  stations  reporting  arrival  times  to  the  ISC.  Other  magnitudes  in  Table  1 
include  the  F-wave  magnitude  for  many  Balapan  UNEs  reported  by 
Ringdal  et  al  (1992),  which  is  a  maximum  likelihood  value  derived  by  the 
British  Atomic  Weapons  Establishment  (AWE);  a  /*-wave  magnitude  also 
obtained  by  AWE  for  many  Degelen  and  Murzhik  UNEs  (personal 
communication  to  Richards  -  these  magnitudes  have  four  significant 
figures)  but  which  is  not  maximum  likelihood;  and  an  Lg-wave  magnitude 
for  many  Balapan  events  reported  by  Ringdal  et  al  (1992).  Finally,  the 
Lg-wave  magnitude  for  two  Balapan  events  (1980  April  25,  and  the  large 
event  of  1980  September  14)  is  given  by  Richards  and  Shi  (1994),  who 
digitized  high  quality  Lg  recordings  from  the  station  WMQ  and  calibrated 
RMS  Lg  measurements  at  WMQ  in  terms  of  the  NORSAR  Lg  magnitude 
scale. 

Note  in  Table  1  that  there  are  ten  events  (reported  by  Bocharov  et  al, 
1989)  for  which  no  magnitude  information  is  listed.  Thus,  the  Table  gives 
locations  for  282  UNEs  (100  at  Balapan,  157  at  Degelen,  25  at  Murzhik), 
and  magnitude  information  for  only  272. 

In  recent  years,  several  Russian  and  Kazakh  sources  have  published 
information  on  the  total  number  of  UNEs  conducted  at  Semipalatinsk  from 
the  first  event  on  1961  October  11,  to  the  last  on  1989  October  19  prior  to 
the  break  up  of  the  U.S.S.R.  and  the  closing  down  of  this  test  site.  For 
example,  Mikhailov  et  al  (1992)  state  that  343  UNEs  were  conducted  at  the 
Semipalatinsk  Test  Site.  And  van  der  Vink  et  al  (1992,  their  Figure  6)  give 
the  number  of  UNEs  in  each  year  at  this  test  site,  the  total  coming  to  343. 

Figure  1  shows  a  histogram  comparing  the  number  of  events  per 
year  from  Table  1,  with  the  number  indicated  by  van  der  Vink  et  al.  Such 
a  comparison  should  be  associated  with  several  caveats,  since  it  is  not 
always  clear  what  constitutes  a  "nuclear  test"  for  counting  purposes.  For 
example,  if  three  nuclear  devices  are  shot  at  the  same  time  in  the  same 
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shaft,  would  that  count  as  one  or  three  tests?  (It  would  be  counted  as  one 
test  according  to  the  definition  of  a  nuclear  test  under  the  verification 
protocol  to  the  Threshold  Test  Ban  Treaty,  and  there  would  be  only  one 
test  apparent  from  seismic  signals.)  If  a  nuclear  test  explosion  was 
initiated  but  did  not  result  in  significant  nuclear  yield,  then  there  would 
typically  be  no  seismic  signal  -  yet  for  some  purposes  this  would  be 
counted  as  a  nuclear  test.  (A  successful  one-point  safety  test  would  have 
these  latter  characteristics,  as  would  a  "fizzle"  —  a  device  that  failed  to 
attain  even  a  small  fraction  of  its  design  yield.  Several  of  the  UNEs 
conducted  in  the  early  years  of  testing  at  the  Nevada  Test  Site  have 
announced  yields  of  either  zero  or  less  than  a  ton  of  TNT  equivalent.)  Note 
in  Figure  1  for  1984  that  there  is  one  more  UNE  in  the  total  of  claimed 
explosions  (based  on  Table  1),  than  in  the  total  given  by  Russian  sources: 
perhaps  one  of  the  claimed  events  is  a  chemical  explosion.  Finally,  the 
number  343  for  Semipalatinsk  may  refer  to  nuclear  weapons  tests  only, 
when  at  least  some  of  the  explosions  at  this  test  site  have  been  described 
in  Soviet  literature  as  part  of  the  Peaceful  Nuclear  Explosion  program. 
Nevertheless,  the  number  reported  from  teleseismic  observations  (the  282 
reported  in  Table  1)  appears  to  be  significantly  lower  that  the  343  now 
reported  as  the  total  carried  out  at  Semipalatinsk,  and  this  difference  was 
a  stimulus  for  our  project  reported  here. 

The  location  of  the  test  site  itself  appears  not  to  have  been  shown 
accurately  in  any  publication  available  in  the  western  literature.  But  when 
the  Threshold  Test  Ban  Treaty  entered  into  force  in  December  1990,  the 
Soviet  Union  and  the  United  States  formally  exchanged  information  on  the 
boundaries  of  their  respective  test  sites.  The  extent  of  the  Semipalatinsk 
Test  Site  was  then  specified  (a)  by  a  list,  given  to  the  U.S.,  of  the  locations 
of  152  posts  on  the  test  site  boundary,  and  (b)  by  a  map,  now  held  in  the 
U.S.  State  Department  Treaty  Library.  There  are  some  slight  discrepancies 
between  (a)  and  (b),  but  the  map  (b)  is  presumably  definitive  and  was 
used  to  prepare  the  illustration  of  test  site  boundaries  shown  in  Figure  2. 

In  order  to  begin  a  systematic  program  of  data  acquisition  for  other 
seismic  events  that  might  or  might  not  be  UNEs  at  Semipalatinsk,  or  that 
were  relevant  to  a  program  of  research  into  regional  discriminants  for  the 
area,  we  first  made  a  list  of  all  the  seismic  events  for  which  we  would  seek 
additional  information.  The  list  appears  here  as  Table  2,  giving  dates, 
sources  of  information,  and  approximate  times,  locations,  and  magnitudes 
(if  available)  for  65  events.  It  includes  the  10  events  in  Table  1  that  lack 
magnitude  information.  It  includes  24  events  pointed  out  by  Sykes  and 
Ruggi  (1989)  and  claimed  by  these  authors  to  be  UNEs,  that  are  located 
only  approximately,  since  often  the  Large  Aperture  Seismic  Array  (LASA) 
in  Montana  was  their  only  source  of  teleseismic  data  for  small  events. 


5 


(Their  list  also  includes  most  of  the  events  in  Table  1.)  The  sources  of 
information  for  Table  2  are  as  follows: 

Boch.  =  Bocharov  et  al  (1989),  10  UNEs  with  no  reported  magnitude; 

H.R.R.  =  Hansen,  Ringdal  and  Richards  (1990),  who  incorrectly  listed 
the  event  of  1988  September  26  as  a  UNE; 

Ring.  =  personal  communication  from  Erode  Ringdal,  who  sent  us 
event  magnitudes  and  estimated  origin  times  for  several 
events  as  detected  by  NORSAR  or  NORESS; 

Ring. (90)  =  Ringdar  (1990),  who  describes  teleseismic  data  for  the 
small  announced  UNE  of  1988  December  28; 

Suit.  =  Djamil  Sultanov  (personal  communication),  who  informed  us 
that  the  event  of  1988  September  26  was  not  a  UNE; 

S&R  =  Sykes  and  Ruggi  (1989); 

V.An  =  Vadim  An  (personal  communication),  whose  information  was 
usually  based  upon  data  recorded  at  the  Borovoye 
Geophysical  Observatory; 

V.Kh.  =  The  working  notes  of  Vitaly  Khalturin,  based  on  observations 
made  in  Central  Asia  prior  to  working  on  this  paper. 


STATIONS  AND  INSTRUMENTATION 

Our  procedure  was  to  study  the  65  events  of  Table  2  using  the 
seismological  data  from  local  networks  receiving  regional  seismic  waves 
from  the  Semipalatinsk  Test  Site.  We  wanted  to  determine  event  locations 
and  magnitudes,  and  to  get  some  evidence  about  their  nature  (whether 
UNE,  earthquake,  or  chemical  explosion). 

We  used  readings  of  records  from  stations  installed  by  the  Complex 
Seismological  Expedition  of  the  Institute  of  Physics  of  the  Earth,  of  the 
Russian  Academy  of  Sciences.  Part  of  this  data  was  obtained  from 
permanent  stations  GRM,  NSB,  TLG,  and  ZRN.  Most  of  the  Expedition's 
stations  were  temporary,  each  working  for  about  2  years.  Also  we  have 
used  the  arrival  times  from  bulletins  of  permanent  stations  in  Central  Asia, 
Kazakhstan  and  Altai,  and  in  some  cases  from  the  digital  station  at 
Borovoye  (station  BRV:  see  Richards  et  al,  1992).  Table  3  lists  the  52 
stations  we  used,  and  Figure  3  shows  a  station  map. 

The  basic  characteristics  of  the  five  different  types  of  instruments 
that  we  used  were  as  follows: 
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Type 

T  —  seconds 

Magnification 

SKD 

0.2  —  18 

(1  —  1.5),  K 

SKM 

0.1  —  1.4 

(30  —  60),  K 

RVZT 

0.2  —  1.2 

(0.1  —0.3),  M 

KSE 

0.7  —  1.1 

(0.5  —  1.0),  M 

ChISS 

multichannel 

0.022  —  40 

and  T2  here  are  the  periods  at  which  the  response  is  down  to  half  the 
power  of  the  peak  response.  Most  stations  were  equipped  with  the 
standard  short-period  instrument  known  as  SKM,  which  records 
displacement  in  a  quite  wide  band  (0.7-10  Hz).  Their  records  or  the 
bulletins  were  used  to  calculate  epicentral  distances  and  origin  time  as 
discussed  below.  Some  of  the  stations  listed  in  Table  3  also  have  the  broad 
band  instrument  SKD,  and  some  had  the  RVZT  instrument.  A  few  stations 
had  the  very  high  gain  narrow-band  KSE  instrument.  All  four  of  these 
seismometers  consist  of  a  coupled  inertial  sensor/galvanometer 
sensor/velocity  pickup,  and  can  be  described  by  the  standard  equations 
for  such  systems  (see  for  example  Aki  and  Richards,  1980,  equation  10.63; 
or,  equation  10.64  since  coupling  is  negligible).  Their  displacement 
responses  are  shown  in  Figure  4a.  It  may  be  a  surprise  to  some  western 
seismologists  that  it  is  possible  to  achieve  such  a  flat  response  to 
displacement  (e.g.  for  SKM  and  SKD)  in  an  instrument  with  a  coil  that 
provides  a  voltage  proportional  to  velocity.  The  way  the  desired  result  is 
achieved,  for  the  SKD  instrument,  is  to  have  a  strongly  overdamped 
galvanometer  (e.g.  Tg^j  =  1  s,  Cg^j  in  the  range  4  to  6)  and  an 
underdamped  inertial  sensor  (e.g.  =  20  s,  £g  in  the  range  0.4  to  0.45), 

with  only  a  small  amount  of  coupling  (less  than  0.1). 

The  fifth  instrument  we  used  is  known  by  its  Russian  acronym, 
ChISS,  and  such  data  were  obtained  from  4  stations:  TLG,  NSB,  ZRN  and 
GRM.  ChISS  was  designed  by  Zapolskii  (1960,  1971)  and  is  described  by 
Zapolskii  and  Khalturin  (1960),  Rautian  and  Khalturin  (1978),  and  Rautian 
et  al  (1978).  The  vertical  component  was  analyzed.  In  this  instrument, 
the  signal  from  the  seismometer,  proportional  to  the  velocity  of  ground 
motion,  is  passed  through  a  system  of  narrow  bandpass  filters  and 
recorded  on  photographic  paper  or  by  ink-recorder.  Thus  the  amplitudes 
on  each  record  are  proportional  to  the  velocity  of  ground  motion  in  a 
particular  frequency  band.  ChISS  system  at  GRM  and  TLG  stations  have  as 
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many  as  16  channels.  We  used  data  mostly  from  8  channels  in  the  range 
from  0.3  to  40  Hz.  For  a  large  event,  (with  magnitude  5  and  more)  the 
Rayleigh  wave  is  strong  enough  to  see,  and  in  that  case  we  widened  the 
frequency  range  studied  to  0.05  Hz  to  include  long  periods.  The  bandwidth 
of  ChISS  channels  Af  is  proportional  to  their  central  frequencies,  /^: 

=  ^/c 

The  parameter  k  is  equal  0.7  in  the  long-period  part  of  the 
frequency  range,  0.48  in  the  central  part  from  0.3  to  10  Hz,  and  0.22  for 
high-frequency  channels,  where  values  are  18,  27  and  40  Hz.  Figure  4b 

(upper)  shows  the  ChISS  spectral  amplitude  response  for  an  8-channel 
system.  The  roll-off  on  either  side  of  each  peak  is  very  steep.  Figure  4b 
(lower)  shows  the  impulse  response  for  four  channels.  Note  that  time  runs 
from  right  to  left,  a  common  convention  in  Russian  seismograms. 

Each  channel  was  calibrated  daily  by  driving  the  pendulum  at 
constant  velocity  amplitude  and  very  slowly  decreasing  frequency.  The 
response  of  this  calibration  signal  recorded  directly  on  the  seismogram  has 
an  envelope  identical  to  the  amplitude-frequency  response  curve.  The 
maximum  amplitudes  of  response  usually  correspond  to  1  |Llm/s  (i.e.  one 
micron/second)  of  velocity  of  ground  oscillation.  This  system  lets  us  get 
the  frequency-dependent  content  as  well  as  the  time-dependent  content  of 
each  seismic  wave.  In  kinematic  studies  we  can  choose  the  channel  where 
the  wave  arrival  looks  most  sharp  and  clear,  to  measure  the  time  arrival 
more  accurately.  Figure  5a  is  an  example  of  a  ChISS  record  showing  six 
channels  (each  with  its  calibration  signal  -  in  this  case  with  maximum 
amplitude  corresponding  to  0.5  pm/s)  with  data  for  the  UNE  of  1976 
December  30  at  Degelen.  Again,  note  that  time  runs  from  right  to  left  in 
each  channel.  This  example  of  ChISS  data  recorded  at  a  distance  of  730 
km  shows  several  different  regional  phases,  as  we  discuss  in  the  section 
below.  Figure  5b  shows  six  channels  of  a  Talgar  ChISS  record  of  the 
Degelen  explosion  of  April  21,  1976.  Five  minutes  after  this  explosion, 
there  was  a  Balapan  explosion  of  similar  magnitude,  for  which  the  Talgar 
ChISS  record  is  shown  in  Figure  5c.  Detailed  examination  of  Figs  5b  and  5c 
shows  that  different  regional  phases  recorded  at  Talgar  are  excited  to  a 
slightly  different  degree  by  Degelen  and  Balapan  explosions. 


REGIONAL  PHASES  IN  CENTRAL  ASIA  AND  KAZAKHSTAN 

Several  different  regional  phases  are  seen  at  distances  from  0  to 
1400  km  from  the  Semipalatinsk  Test  Site. 
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At  distances  up  to  about  240  km  the  main  regional  phases  are  P  and 
s.  Their  velocities  are  6.2  and  3.54  km/s  respectively.  Beyond  230-250 
km  the  Pn  wave  appears  as  a  first  arrival.  This  wave  has  stronger  high- 
frequency  spectral  content  than  all  other  regional  waves. 

As  shown  by  Nersesov  and  Rautian  (1964),  the  Pn  wave  is  not 
simple.  It  consists  of  2  or  3  waves.  Each  successive  wave  appears  as  a 
later  arrival,  with  amplitude  larger  than  the  previous  one.  As  distance 
increases,  the  Pn2  time  arrival  moves  closer  to  the  first  arrival  and 
amplitudes  of  Pn2  decrease.  We  do  have  not  a  good  range  of  distances 
with  the  data  of  this  paper  to  see  this  picture  clearly.  We  can  say  only  that 
the  interval  of  distances  where  we  observe  Pn2  is  small  -  between  500 
and  750  km.  Pn2  is  observed  at  some  stations,  but  not  at  each  of  them. 

At  TLG  the  Pn2-wave  arrives  about  8-10  s  later  then  Pni  and  has  a 
larger  amplitudes  (2-5  times).  Examples  can  be  seen  in  Figures  5abc  (the 
phases  marked  Fj  and  ^2)-  ZRN  the  amplitudes  of  Pn2  are  small.  At 
NSB  (all  three  stations  are  near  the  same  distance  from  the  test  site),  Pn2 
is  absent.  The  Pn^  amplitudes  attenuate  up  to  1600-1800  km,  and  then 
increase  with  distance,  presumably  due  to  a  deep  mantle  arrival  (which, 
however,  is  a  continuation  of  the  straight-line  travel-time  curve,  with 

only  a  slight  increase  in  phase  velocity).  The  amplitude  behavior  is 
strongly  different  depending  on  the  direction  of  wave  propagation  from 
the  Semipalatinsk  Test  Site.  At  the  East  direction  amplitudes  decrease 
with  distance  steeper  and  then,  after  1600-1800  km,  increase  sharper, 
than  in  the  West  direction. 

The  Pg  wave  can  be  observed  up  to  distances  of  about  700-800  km. 
The  Pg  arrival  is  impulsive  at  distances  250-400  km.  Its  amplitudes  are 
more  than  Pn  (5-10  times).  But  beyond  500  km  Pg  looks  like  a  group  of 
low-frequency  waves  with  an  emergent  arrival.  Its  amplitudes  are  of  the 
same  order  as  P-coda  amplitudes,  increasing  gradually  with  lapse  time. 

The  Sn  wave  diverges  from  Lg  near  280  km.  It  can  be  seen  clearly 
only  on  40-50%  of  the  records.  If  the  propagation  path  is  through  the 
shield,  Sn  is  stronger  and  can  be  seen  better  then  in  an  orogenic  area.  In 
the  shield  it  has  nearly  the  same  high-frequency  spectral  content  as  Pn. 

Beyond  1200-1400  km  the  S  wave  appears  about  20-25  s  later  then 
the  Sn  wave  (unlike  the  relation  between  Pn  and  P,  for  which  there  is 
no  time  offset).  This  is  the  S  wave,  with  phase  velocity  5.7  km/s  pointed 
out  in  J-B  tables.  It  arrives  clearly  after  1400-1500  km.  Its  spectral 
content  is  much  more  low-frequency,  than  that  of  Sn.  At  Garm  (1350  km 
from  the  test  site)  one  can  observe  both  S  and  Sn  on  the  records  from 
Semipalatinsk.  The  regional  variations  in  travel  time  of  Sn  and  especially 
for  S  at  1300-1500  km  are  very  large:  up  to  about  8  s. 

Before  the  Lg  wave  arrival  we  sometimes  find  the  low-frequency 
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wave  called  Li. 

At  all  the  distances  under  consideration  the  dominant  wave  is  Lg.  It 
appears  to  consist  of  several  groups,  from  which  we  can  identify  arrivals 
that  we  call  Lg^  and  Lg2-  The  Lg^  arrival  is  impulsive,  whereas  Lg2  is 
emergent,  with  increasing  amplitudes  within  a  long  group  of  interfering 
Lg  oscillations. 

The  final  recognizable  wave  is  a  Rayleigh  wave,  Rg,  but  if  the 
magnitude  of  the  explosion  is  less  then  4,  the  signal-to-noise  ratio  is  too 
small  for  this  wave  to  be  seen.  This  point  requires  some  discussion.  Thus, 
the  terminology  here  implies  a  dispersive  wave,  but  in  practice  there  is 
little  surface  wave  dispersion  over  paths  extending  less  than  500  km  from 
Semipalatinsk:  periods  from  about  5-6  s,  up  to  25-30  s,  are  superimposed. 
The  maximum  amplitude  of  Rg  for  an  explosion  recording  appears  to  have 
a  period  of  about  5  to  6  s;  for  an  earthquake,  the  period  is  longer.  Note 
that  for  a  conventional  measurement  of  surface-wave  magnitude,  M 
made  from  a  fundamental  mode  20- s  Rayleigh  wave  at  a  distance  where 
the  wave  is  well  dispersed,  the  difference  between  body-wave  magnitude 
nify  and  surface  wave  magnitude  for  a  nuclear  explosion  is  on  average 
about  1.2.  But  for  a  measurement  made  around  500  km,  where  M ^  is 
measured  at  5  to  6  s  period,  the  difference  is  about  0.6. 

However,  it  needs  a  fairly  large  explosion  (say,  >  4.8)  for  such  Rg 
observations  to  be  made. 

Although  we  had  access  to  some  data  from  over  50  stations,  in 
practice  our  data  came  mostly  from  about  ten  stations.  At  the  other 
extreme,  about  60%  of  the  stations  were  used  for  only  1-3  events.  The 
number  of  events  for  which  each  station  was  used  was  as  follows: 


# 

Sta. 

# 

Sta. 

# 

Sta. 

# 

Sta, 

# 

Sta. 

# 

Sta. 

# 

TLG 

60 

KRG 

9 

4 

BOM 

2 

SVE 

2 

NGN 

B 

FRZ 

7 

4 

CHGU 

2 

TDK 

2 

DZK 

1 

PDG 

1 

UKN 

31 

KRR 

6 

4 

CHK 

2 

TGM 

2 

B 

RYB 

1 

BRV 

25 

MDO 

6 

ATA 

3 

KBZ 

2 

B 

TUR 

1 

NSB 

23 

PRZ 

6 

KRD 

3 

KKB 

2 

KKR 

1 

VED 

1 

KRM 

16 

SEM 

6 

3 

MIX 

2 

ALB 

1 

KSU 

1 

ZRN2 

11 

MAK 

5 

ULG 

3 

MRT 

2 

KUS 

1 

BYK 

4 

VOS 

3 

SKL 

2 

ART 

1 

MKR 

1 
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The  availability  of  quiet  sites  in  the  region  enables  the  deployment 
of  high  gain  instrumentation.  The  RVZT  and  KSE  records  of  a  presumed 
UNE  at  Degelen  on  January  29,  1971,  are  shown  in  Figure  6ab,  indicating 
excellent  signal- to-noise,  in  this  example  with  gains  of  120,000  and 
1,100,000. 


THE  TRAVEL-TIMES  IN  A  REGION  AROUND  THE  TEST  SITE 

Most  of  the  events  we  have  studied  are  small  and  high-gain  records 
were  not  often  available.  Thus,  in  some  records  the  Pn^  wave  cannot  be 
found  and  the  only  arrival  times  that  can  be  read  are  for  shear  waves, 
such  as  Sn,  Li  and/or  Lg.  So  we  need  to  use  such  shear  wave  data  when 
calculating  the  epicentral  distances  and  origin  time.  Earthquakes  in  the 
nearby  orogenic  Tien  Shan  region  were  studied  by  Nersesov  and  Rautian 
(1964)  and  Shazilov  (1989),  but  in  our  case  the  seismic  waves  propagated 
in  the  Kazakh  platform  and  we  needed  to  obtain  the  travel  times  versus 
distances  for  this  zone.  For  this  purpose  we  used  the  data  from  UNEs  for 
which  the  epicenters  are  known  very  accurately. 

The  Pn^  velocity  is  practically  constant  up  to  1200-1300  km  and  is 
8.1  km/s,  and  then  slightly  increasing  up  to  8.5  km/s  over  the  range 
1300-1700  km. 

The  velocity  of  Pn2  is  8.3  km/s. 

The  velocity  for  P  as  well  as  for  Pg  is  the  same,  6.2  km/s.  This 
value  is  a  little  higher  than  usually  observed  in  the  Tien  Shan  and  Pamir 
regions,  and  stays  constant  throughout  the  distance  range  within  which 
these  waves  are  observed. 

Figure  7  shows  the  travel  times  for  9  different  regional  phases  in 
East  Kazakhstan.  The  following  relationships  between  distance  and  arrival 
time  are  available  to  determine  epicentral  distances  D  in  kilometers: 


D  =  6.2*  t(p)_ 

D  =  3.54  *  t(s) 

D  =  8.1  *  t{Pn{)-69 
D  =  8.3  *  t{Pn2)-\66 
D  =  6.1*  t(Pg) 

D  =  4.62  *  t(Sn)-69 
D  =  3.9  *  tiLi) 

D  =  3.54  *  tiLgi) 

D  =  3.4*  t(Lg2) 


0  -  240  km 
0  -  240 
240  -  1200 
500  -  750 
240  -  900 
240  -  1400 
600  -  1000 
240  -  1400 
600  -  1400. 
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The  reliability  of  measured  arrival  times  for  small  events  is  poor  due 
to  the  small  signal-to-noise  ratios,  and  results  in  a  typical  standard  error  of 
about  15  km  when  the  above  relations  are  used  to  estimate  D.  The  error 
in  estimated  epicentral  distances  is  smallest  for  Pn^  and  Lg^,  when  D 
calculated  from  the  above  equations  is  compared  with  the  real  distances 
from  published  accurate  epicenters. 


DETERMINATION  OF  EPICENTERS 

By  taking  the  difference  in  arrival  time  between  pairs  of  phases  at  a 
single  station,  and  using  the  above  relationships  between  distance  and 
arrival  time,  we  can  estimate  the  epicentral  distance  from  each  station  that 
has  adequate  data.  The  actual  epicenter  can  then  be  found  graphically  as 
the  point  of  intersection  of  circles  with  the  stations  at  their  centers  and 
radius  being  the  epicentral  distance  from  the  station.  Figure  8  shows 
several  examples.  Using  the  data  from  many  stations,  we  get  an 
intersection  point  from  each  pair  of  stations.  The  epicenter  is  calculated  by 
averaging  all  points  of  intersection,  with  a  weight  for  each  point,  equal  to 
the  sine  of  the  intersection  angle. 

If  the  epicenter  is  close  to  the  center  of  the  Test  Site,  we  can  simplify 
the  calculation  of  epicenters,  using  straight  lines  rather  than  arcs  of  circles. 
To  do  this  we  choose  the  central  point  of  Degelen,  (49.85°E,  78.08°N),  as  a 
master  epicenter  and  estimate  the  value  of  AD  =  D  -  Dq  where  Dq  is 
the  distance  from  the  station  to  the  central  point.  In  many  cases  AD  is 
only  a  few  percent  of  Dq.  For  example,  for  the  event  of  1988  December  28 
we  had  several  phases  recorded  at  7  different  stations.  The  average  value 
of  AD,  and  the  value  of  Dq,  were  as  follows  (all  in  km): 


AD 

^0 

E 

KSU 

-4.0 

21S 

E 

SEM 

4.8 
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E 

UKN 

-6.4 

488 

S 

MRT 

9.1 

525 

S 

PDG 

-9.6 

727 

s 

TLJ3 

4.2 

733 

w 

m£ 

0.7 

211. 

Letters  E,  N,  S  or  W  here  indicate  the  direction  from  source  to  station 
(whether  East,  North,  South  or  West).  Because  AD  is  so  much  smaller  than 
Dq,  the  curvature  of  arcs  is  low  and  we  can  replace  the  arcs  by  straight 
lines,  perpendicular  to  the  direction  from  epicenter  to  stations. 
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If  data  were  available  from  only  2  stations,  we  can  get  a  single 
intersection  of  straight  lines.  The  longer  axis  of  the  error  area  was 
calculated  as  10  km  divided  by  sin(azl-az2),  and  the  minimum  was  taken 
as  10  km.  The  azimuth  of  max  error  was  taken  as  the  average  azimuth  of 
the  pair  of  stations. 

Several  UNEs  with  accurate  epicenters  were  used  to  check  the 
accuracy  of  epicenter  determination.  Results  are  shown  in  Table  4. 
Calculations  were  done  in  2  versions:  using  only  Pn^  (indicated  as  Pn  in 
the  Table),  and  using  the  average  error  in  location  from  all  waves  where 
more  than  just  Pn  was  used.  In  most  cases  the  error  is  about  1-12  km, 
but  for  a  few  UNEs  it  is  large,  up  to  100  km.  The  low  accuracy  can  be  due 
to  error  of  arrival  times  (for  some  cases  only  bulletin  data  were  available), 
and/or  because  stations  were  available  only  in  one  general  direction  (e.g. 
South)  from  the  source.  The  standard  deviation  for  an  estimate  of  the 
distance,  using  Pn^  for  a  single  station,  is  12  km. 

One  can  see  from  Table  4  that  the  accuracy  of  epicenter  location  is  of 
the  same  order  for  Pn  alone,  as  for  averaged  data  from  all  available 
waves.  So,  if  Pn  data  of  many  stations  are  available,  we  need  not  use 
other  waves.  But  for  small  events,  when  Pn  is  not  readable  on 
seismograms,  we  used  all  wave  data  available. 

The  list  of  52  stations  (see  Table  3)  includes  10  to  the  East,  1  to  the 
North,  24  to  the  South,  and  17  to  the  West.  If  the  stations  recording  an 
event  are  all  from  only  one  side,  (for  example  May  7  and  Aug  19,  1966, 
where  we  had  data  only  from  the  South  side),  then  only  one  coordinate  can 
be  found  accurately.  The  number  of  stations  for  which  data  are  available 
varied  for  each  event  from  1  (11%),  2-4  (46  %);  5-9  (27  %)  ;  to  10-11  (only 
8  %).  For  10  events  the  data  were  obtained  from  stations  localized  only  to 
one  side,  usually  South  or  East. 


DETERMINATION  OF  EPICENTERS 

Of  the  65  events  studied  in  this  paper  (Table  2),  10  already  have 
known  locations  (Bocharov  et  al,  1989).  In  this  section  we  describe  our 
first  location  estimates  for  the  remaining  55  events,  presenting  the  results 
in  Tables  5,  6,  7,  and  8  (and  see  also  Figure  8,  for  12  events).  Throughout 
the  year  in  which  we  worked  on  this  project  together  at  Lament,  we 
continued  to  receive  additional  information  on  arrival  times  for  specific 
events  at  specific  stations.  Table  12,  in  a  later  section  presenting  our  final 
results,  reflects  this  additional  information,  and  reports  hypocentral 
estimates  that  are  slightly  different  for  some  events  from  the  results 
reported  in  this  section. 
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Table  5  gives  our  first  location  estimates,  for  events  that  appear  to 
be  on  or  very  close  to  the  Semipalatinsk  Test  Site,  and  for  which  we  have 
data  available  from  more  than  one  general  azimuthal  direction.  These  are 
our  best  located  events. 

Table  6  describes  four  events  for  which  data  were  available  at  only 
one  azimuth,  either  to  the  South  (three  events)  or  the  East  (one  event).  In 
these  cases  we  can  estimate  only  one  coordinate,  but  it  is  consistent  with  a 
location  on  the  test  site.  The  estimated  origin  time  in  three  of  these  cases 
is  consistent  with  being  on  the  exact  minute  (00  seconds)  -  further 
evidence  that  these  events  are  explosions. 

Table  7  describes  six  events  for  which  data  were  available,  at  only 
one  station.  The  calculated  epicentral  distances  D  for  two  of  the  events 
are  practically  the  same  as  the  distance  Dq  io  the  central  point  of  Degelen. 
For  four  of  the  events,  AD  =  D  -  Dq  was  as  much  as  10-30  km. 

Although  we  cannot  find  these  epicenters,  the  data  do  not  contradict  our 
assumption  that  these  events  lie  on  the  test  site. 

Table  8  describes  three  events  (1967  June  3,  1967  July  16,  1981 
March  31)  at  distances  100-300  km  from  the  test  site;  and  three  more 
events  that  are  even  farther  away.  Thus,  the  event  of  1980  October  26  is 
1545  km  from  Talgar,  more  than  twice  the  distance  from  Talgar  to  Degelen 
(733  km).  The  event  of  1983  February  13  is  765  km  from  station  UKN, 
whereas  the  distances  from  UKN  to  Degelen  and  Balapan  are  488  and  425 
km  respectively  The  event  of  1984  August  26,  for  which  we  have 
readings  from  two  stations,  lies  at  one  of  the  two  intersection  of  two 
circles.  Both  choices  of  epicenter  lie  far  from  the  test  site. 

Commenting  in  general  on  this  work  of  event  location,  we  can  say 
that  the  network  of  regional  stations  in  Central  Asia  uncovered  many 
events  additional  to  those  located  from  teleseismic  data.  The  most 
important  stations  are  those  at  distances  less  that  500  km  from  the 
Semipalatinsk  Test  Site.  Even  at  700  km,  detection  becomes  significantly 
worse.  The  worst  distances  are  around  1000  km.  At  2000  km  and 
beyond,  amplitudes  are  stronger  than  at  1000  km,  these  signals  coming 
from  the  deep  mantle  phases  P  and  S  ,  and  teleseismic  methods  apply. 
These  results  about  amplitude  dependence  on  distance  have  been  known 
for  more  than  30  years,  and  were  well  understood  during  test  ban 
negotiations  in  the  period  1960-1963.  The  data  of  Sykes  and  Ruggi  (1989) 
for  small  events,  based  principally  upon  LASA  signals,  indicated  that 
teleseismic  data  can  be  useful  for  detection  of  events  with  magnitude 
around  4  and  below,  but  that  locations  based  upon  data  from  a  single 
teleseismic  array  are  often  poor  (compare  Sykes'  and  Ruggi's  locations 
noted  in  Table  2,  with  our  locations  in  Tables  5-8). 
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DETEilMIN ATION  OF  MAGNITUDES 

We  found  there  are  significant  differences,  some  of  which  are 
systematic,  between  the  determinations  of  magnitude  by  ISC,  NORSAR,  and 
AWE  (the  British  Atomic  Weapons  Establishment)  for  small  events.  The 
ISC  and  AWE  magnitudes  are  both  based  upon  teleseismic  F-wave 
observations;  and  the  NORSAR  magnitude  is  based  upon  RMS  Lg 
measurements  (which  usually  require  Semipalatinsk  events  to  be  above 
magnitude  5.5  in  order  to  have  an  adequate  signal- to-noise  ratio  at 
NORSAR,  more  than  4000  km  away).  It  seems  to  us  that  the  AWE 
determinations  are  preferable,  especially  for  events  with  magnitudes  less 
then  5.  The  differences  between  AWE  and  ISC  are  due  to  the  systematic 
error  that  can  follow  when  working  with  signals  that  are  close  to  the  noise 
level,  since,  in  the  ISC  procedure,  signals  that  are  below  the  noise  level  are 
not  allowed  for  in  forming  the  average,  so  the  averaged  magnitude  is 
biased  high.  The  AWE  maximum  likelihood  magnitude  allows  for  the 
statistics  of  non-observation  in  the  presence  of  noise.  Figure  9a  shows  a 
comparison  of  ISC  and  AWE  magnitudes,  and  the  ISC  magnitudes  are 
indeed  to  be  too  high  for  the  smaller  events.  Figure  9b  shows  that  the 
AWE  and  NORSAR  magnitudes  are  in  good  agreement,  as  noted  and 
discussed  by  Ringdal  et  al  (1992). 

However,  although  we  prefer  AWE  values,  we  recognize  that  ISC 
magnitudes  are  widely  used  as  a  standard.  Therefore  we  have  used  ISC 
mf^(P)  and  separately  AWE  and  NORSAR  magnitudes  when  comparing  our 
regional  measurements  with  teleseismic  observations. 

We  used  two  ways  to  calculate  magnitude  from  regional  data.  One 
way  is  based  on  the  K  scale  (energy  class)  which  is  used  in  the  former 
Soviet  Union  in  all  local  networks.  K  is  approximately  equal  to  log  E 
where  E  is  the  radiated  energy  in  joules  (Rautian,  1960).  To  derive  the 
relationship  between  K  and  teleseismic  magnitude  and  to  calculate  the 
magnitudes  Af(K)  reported  in  Table  9,  we  took  K  values  from  Soviet 
bulletins  or  calculated  K  ourselves  from  amplitude  readings.  Using  the 
data  on  K  and  the  ISC  values  of  /n^(F)  for  large  UNEs  at  the  Semipalatinsk 
Test  Site  the  relationship  can  be  found  as 

m^(P)  =  0.42  K. 

The  values  of  ISC  m^(F)  are  not  good  for  small  events  (as  noted  from 
Figure  9a),  and  if  we  compare  K  with  the  magnitudes  from  AWE  or  NORSAR 
we  find  a  different  relationship: 

m(AWE  or  NORSAR)  =  0.44  K  -  0.53. 
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The  points  to  which  this  last  equation  is  fit,  are  shown  in  Figure  9c  across 
six  orders  of  magnitude  in  energy.  (In  addition  to  the  m^(P)  values  from 
AWE  and  the  m^(Lg)  values  from  NORSAR,  we  sometimes  used  m^(P) 
values  from  NORSAP.  for  small  events.)  Either  of  the  above  two  formulae 
may  be  used  to  estimate  a  magnitude  based  on  K,  denoted  Af(K).  We 
preferred  to  use  the  second  formula. 

Our  next  method  of  determining  m^(P)  was  to  get  the  separate 
relationship  between  log  A  (from  ChISS  records  at  TLG,  NSB  and  ZRN)  and 
ISC  mfjiP)  for  each  station  using  the  large  events.  We  then  extrapolated 
the  formula  using  the  observed  values  of  log  A  to  estimate  magnitudes  for 
the  smaller  events. 

The  equations  were  obtained  separately  for  each  wave  (P«i,Pn2, 

Pg,  Sn,  Lg^,  Lg2,  Rg  and  amplitudes  of  coda  at  the  lapse  time  500 
seconds).  Also  we  used  the  averaged  amplitudes  of  all  P  waves;  the  sum 
of  amplitudes  of  Lgj,  Lg2,  and  Lg^,',  and  the  absolute  max  of  Lg, 
independent  of  its  time  of  arrival.  The  coda  is  the  best  parameter  for 
magnitude  calculation,  giving  the  more  accurate  results.  But  for  the  small 
events  we  cannot  see  coda  on  records  late  enough  in  lapse  time. 

The  estimations  were  made  for  each  frequency  band  separately.  In 
general  the  frequency  range  of  ChISS  used  in  this  problem  was  from  14  s 
to  2.5  Hz.  The  correlation  equations  obtained  are  as  follows: 

m^(P)  =  /:*log  A  +  b. 

with  the  values  of  k  and  b  being  different  for  each  station,  each  wave, 
and  each  frequency  band.  Figure  9d  shows  the  comparison  between 
m^(AWE)  and  determined  from  ChISS  measurements. 

Small  events  sometimes  were  recorded  by  only  2-3  channels, 
corresponding  to  the  maximum  of  the  spectrum,  and  coda  cannot  be  seen 
and  used.  So  to  estimate  a  magnitude  for  a  small  event  we  can  use  only  a 
limited  number  of  the  above  correlations  between  m^(P)  and  ChISS 
amplitudes,  corresponding  to  those  for  which  we  have  ChISS  data.  We 
form  an  average  over  all  the  magnitudes  derived  from  ChISS  readings  in 
this  way.  The  system  of  coefficients  in  the  last  equation,  above,  to 
calculate  m^(P)  from  ChISS  amplitudes,  was  obtained  using  ISC 
magnitudes.  To  get  a  second  version  of  m^(ChlSS),  corresponding  to  AWE 
and  NORSAR  estimations,  we  use  the  relationship  (e.g.  from  Figure  9a) 

m^(AWE  and  NORSAR)  =  1.25  m^(ISC)  -  1.50  (for  Balapan) 

=  1.25.m^(ISC)  -  1.37  (for  Degelen). 
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Our  magnitude  results  are  given  in  Table  9,  and  we  see  that  values 
for  previously  unknown  events  typically  lie  in  the  range  3.5  —  4.5,  so  these 
explosions  are  mostly  subkiloton.  The  magnitude  of  the  first  Semipalatinsk 
UNE,  assigned  here  as  4.81,  is  somewhat  uncertain.  This  event  (on  October 
11,  1961)  had  an  amplitude  6.7  times  smaller  than  that  of  the  second  UNE 
(on  February  2,  1962)  on  the  Talgar  record  -  and  the  latter  event  had  a 
relatively  well-determined  K  value  of  14.0  from  four  stations  and  hence  a 
magnitude  of  5.63.  Thus  the  magnitude  of  the  first  UNE  at  the 
Semipalatinsk  Test  Site  is  assigned  as  5.63  -  log  6.7  =  4.81. 

In  view  of  the  often  excellent  teleseismic  detection  capability  of 
arrays  (Ringdal,  1990),  it  was  a  surprise  to  us  that  the  events  were  in 
many  cases  not  detected  by  NORSAR.  One  difficulty  with  assessing  the 
performance  of  teleseismic  arrays,  with  our  event  set,  is  that  until  about 
1980  it  was  not  common  for  array  data  to  be  archived  on  a  continuous 
basis.  Only  the  detections  that  were  noticed,  were  saved.  (In  the  case  of 
LASA,  not  even  detections  were  saved,  except  within  a  limited  period  of 
operation  of  this  array.)  We  checked  with  Yellowknife  and  learned  that 
the  only  events  confirmed  at  tliat  array  after  1978  (when  it  is  possible  to 
check  against  an  archive)  were  those  of  1981  March  31  and  1982  June  11 
(personal  communication  from  Bob  North).  This  failure  to  demonstrate 
teleseismic  detections  is  important  to  explore  further,  since  plans  are 
under  development  to  build  a  global  network  of  so-called  alpha  stations, 
mostly  arrays,  intended  to  provide  continuous  and  virtually  complete 
coverage  of  Eurasia  at  the  magnitude  3.5  level  and  better  using  teleseismic 
signals,  for  purposes  of  monitoring  compliance  with  a  Comprehensive  Test 
Ban  Treaty  (CTBT).  If  array  detection  is  not  reliable  at  the  magnitude  3.5 
level  using  teleseismic  signals,  then  increased  attention  must  be  paid  to 
regional  signals  for  purposes  of  detection  as  well  as  identification  -  a 
conclusion  that  has  implications  for  the  numbers  of  seismic  stations  needed 
in  an  effective  network  for  CTBT  monitoring. 


DISCUSSION  AND  CONCLUSIONS 

The  problem  of  identification  of  seismic  events  at  the  Semipalatinsk 
Test  Site  is  considerably  complicated  by  the  occurrence  of  chemical 
explosions.  To  begin  the  discussion  of  identification.  Table  10  reports  for 
this  test  site  the  numbers  of  UNEs  each  year  claimed  by  van  der  Vink  et  al 
(1993),  the  number  of  UNEs  each  year  that  have  been  identified  from 
teleseismic  data  (see  Table  1),  and  the  number  of  chemical  explosions 
greater  than  10  tons  in  each  year  from  1970-1988  reported  by  Saybekov 
(1993).  The  Natural  Resources  Defense  Council  and  the  USSR  Academy  of 
Sciences  set  off  a  10  ton  chemical  explosion  near  the  test  site  in  1987,  that 
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was  recorded  successfully  at  Talgar  (a  distance  of  more  than  700  km),  so 
we  should  expect  no  problem  with  observability  of  chemical  explosions  at 
regional  distances  in  Central  Asia 

One  important  indicator  that  an  event  was  an  explosion,  is  whether  it 
was  shot  on  the  exact  minute.  Prior  to  September  1980,  it  appears  that 
virtually  all  previously  known  UNEs  at  Semipalatinsk  were  carried  out 
either  on  or  very  close  to  an  exact  minute.  The  only  exceptions  in  Table  1 
are  a  single  UNE  at  Murzhik  (1978  July  18),  and  five  UNEs  that  in  each  case 
were  the  second  event  of  a  pair  carried  out  a  few  seconds  apart,  with  the 
first  shot  being  on  the  minute  (1972  December  10,  1974  January  30,  1977 
October  29,  1978  August  29,  and  1978  November  29).  For  many  years,  a 
practical  rule  for  identifying  a  UNE  at  Semipalatinsk  was  to  note  any  event 
whose  first  arrival  was  about  29  s  after  the  minute  at  Borovoye,  and  38  s 
after  the  same  minute  at  Talgar!  We  find  that  many  of  the  previously 
unknown  events  studied  in  this  paper  did  occur,  if  prior  to  September 
1980,  on  or  very  close  to  an  exact  minute.  However,  this  method  does  not 
work  after  that  date,  nor  does  it  help  us  directly  to  distinguish  between 
chemical  and  nuclear  explosions. 

We  have  not  carried  out  a  systematic  study  of  spectral  ratios  (using 
ChISS  data),  to  attempt  to  discriminate  between  earthquakes  and 
explosions,  although  a  preliminary  study  of  such  ratios  has  indicated  that 
P  to  5  amplitude  ratios  are  somewhat  sensitive  to  event  type,  when 
signal-to-noise  ratios  are  high.  We  did  generate  a  preliminary  summary  of 
our  results  for  48  events  -  locations,  magnitudes,  and  a  brief  and 
somewhat  subjective  comment  on  event  type  -  which  is  given  here  as 
Table  11.  Almost  all  of  these  events  are  listed  as  explosions.  The  only 
exceptions  are  the  events  of  1966  December  26  (listed  as  an  earthquake 
because  of  the  time  of  day  and  location  just  off  the  test  site),  1980 
November  6  (listed  as  an  earthquake  in  this  preliminary  summary  because 
of  the  time  of  day),  1988  September  26  (which  a  Russian  colleague  had 
told  us  was  an  earthquake  and  not  a  UNE  -  but  which  we  have  listed  as 
either  an  earthquake  or  a  chemical  explosions),  and  1989  October  20 
(which  a  Russian  colleague  had  told  us  was  a  UNE  collapse). 

Note  in  Table  11  that  prior  to  1978  the  explosion  candidates  are 
listed  as  UNEs  (except  for  one  known  ChE  in  1961),  but  for  1978  and  later 
they  are  listed  as  UNE  or  ChE  (except  for  two  UNEs  at  Degelen).  The  reason 
for  this  change  in  listing  is  that  we  noticed  several  differences  between 
events  occurring  before  and  after  about  1979,  as  follows; 

•  the  error  in  latitude  determination  became  worse  with  time 
(for  1964-1979,  out  of  23  events  only  4  have  latitude  error  >  10  km; 
for  1980-1988,  out  of  18  events,  only  4  have  latitude  error  <  10  km). 
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•  We  developed  our  travel  time  table  and  identification  of  wave  types 
using  the  well-recorded  data  from  known  UNEs  from  the  test  site.  Events 
from  other  locations  can  have  slightly  different  wave  properties.  For  all 
events  up  till  1977  for  which  we  have  the  records,  the  wave  form  is 
approximately  similar.  But  after  1978  at  least  half  of  the  events  have  a 
slightly  different  wave  form:  the  arrival  of  the  Lg  group  is  not  so  sharp; 
and  there  is  a  very  impulsive  arrival  of  Pn2  on  the  ChISS  5  Hz  channel. 

We  called  such  events  the  Balapan  series. 

•  Before  1973  all  our  small  events  occurred  at  Degelen.  In  1973-1978 
only  half  occurred  at  Degelen.  In  1979-1984,  among  19  events  only  one 
could  have  occurred  at  Degelen.  And  the  last  six  events  were  again  mostly 
at  Degelen. 

•  For  the  period  1981-1984  the  number  of  UNEs  claimed  per  year  by 
Russian  sources  (see  Figure  1  or  Table  10)  was  about  equal  to  the  number 
of  events  already  identified  as  UNEs  (Figure  1  or  Table  1),  so  our  small 
events  for  this  period  appear  to  be  in  excess  of  the  number  of  UNEs. 

We  therefore  suspected  at  the  time  of  preparing  Table  11  that  some  of  the 
Balapan  series  were  special  chemical  explosions  -  quite  large  ones,  in  view 
of  their  seismic  magnitudes. 

When  we  circulated  this  preliminary  sununary  table  to  a  small 
number  of  individuals  for  comment,  we  received  the  following  additional 
information  from  Jack  Murphy: 

"Of  the  21  eveats  since  1970  v\Mch  you  locate  as  being  at 
or  near  Balapan  [l^le  11  identifies]  4  as  probable  USE 
(probability  >  0.9) ,  9  as  likely  USE  (0.5  <  probability  < 

0.9) ,  1  as  definite  BQ  (probability  =  1.00) ,  6  as  likely  CE 
(probability  >0.7)  and  1  as  about  equally  likely  to  be 
either  CE  of  BQ.  In  fact,  more  reliable  data  indicate  that 
none  of  these  events  are  USE  or  BQ. " 

Obviously,  this  comment  reinforced  our  opinion  that  a  significant  fraction 
of  the  events  we  have  studied  are  chemical  explosions  at  Balapan.  The 
above  statement  can  be  used  to  infer  that  21  events  in  Table  11,  on  or 
near  Balapan,  are  chemical  explosions.  We  continue  to  believe  that  the 
explosions  we  have  found  at  Degelen  are  UNEs. 

Before  giving  our  final  results.  Figure  10  shows  the  location  of  282 
UNEs  at  Semipalatinsk,  one  earthquake,  and  one  chemical  explosion,  all 
previously  known  (further  information  for  these  UNEs,  given  in  Table  1). 

Our  final  results  are  given  in  Table  12  and  Figures  11  and  12.  Figure 
11  shows  our  locations  for  43  events  not  previously  known  on  the 
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Semipalatinsk  Test  Site.  Most  of  those  at  Degelen,  we  believe  to  be  UNEs. 
Most  of  those  at  Balapan,  we  now  believe  to  be  chemical  explosions.  (Two 
of  the  ChE  events  are  co-located  in  Figure  11,  as  are  two  pairs  of  UNEs  on 
Degelen.)  Figure  1?  shows  events  in  our  study  that  are  located  off  the  test 
site,  but  are  still  within  about  300  km. 

It  should  be  emphasized  that  we  have  no  seismological  basis  for 
discriminating  between  UNEs  and  chemical  explosions.  We  have  chosen  to 
make  the  final  distinction,  in  Table  12,  by  inference  from  Murphy's 
comment  on  our  preliminary  Table  11  (which  named  many  events  at  or 
near  Balapan  as  either  UNE  or  ChE).  Note  that  we  used  seismic  methods  to 
detect  them  -  and  indeed  to  label  them,  in  almost  all  cases  -  as  some  type 
of  explosion.  But  under  a  CTBT,  it  would  appear  to  be  necessary  to  use 
non-seismic  methods  (for  example,  on-site  inspection)  to  confirm  that 
explosions  as  large  as  these  did  indeed  have  a  chemical  nature. 

One  particularly  interesting  event  is  that  of  December  26,  1966, 
which  we  believe  to  be  a  small  earthquake  just  off  the  test  site  to  the 
southeast  (see  Figure  12).  Figure  13a  shows  the  three-component  SKM 
record  of  this  event  at  a  station  only  260  km  away.  The  P  and  S  waves 
both  have  very  high  frequency  content,  and  Lg  has  a  sharp  onset  and  an 
amplitude  several  times  bigger  than  that  of  Pg  (which  in  turn  is  more 
than  20  times  bigger  than  Pn).  Figure  13b  shows  a  similar  record  of  a 
presumed  small  UNE  at  Degelen,  at  a  distance  of  306  km.  The  regional 
phases  are  different  in  several  respects  from  those  of  Figure  13a  -  for 
example,  in  the  Pn  to  Lg  ratio. 

Our  list  of  UNEs  in  Table  12  includes  13  occurring  in  the  period  prior 
to  1973  that  was  the  subject  of  the  paper  by  Bocharov  et  al  (1989).  We 
therefore  believe  that  the  96  UNEs  Bocharov  et  al  describe  were  not  a 
complete  set  of  UNEs  for  1961-1972  at  Semipalatinsk. 

The  problem  remains,  of  a  significant  difference  between  the  303 
UNEs  we  would  now  claim  at  Semipalatinsk  (282  from  Table  1,  plus  21 
additional  UNEs  from  Table  12),  and  the  343  claimed  as  the  total  by 
several  Russian  sources  (see  earlier  discussion:  no  doubt  some  very  small 
UNEs  at  Semipalatinsk,  as  at  the  Nevada  Test  Site,  did  not  generate 
observable  seismic  signals,  even  at  regional  distances). 

In  conclusion,  we  believe  we  have  successfully  located  and  identified 
about  50  small  events  on  or  near  the  Semipalatinsk  Test  Site,  not 
previously  listed  or  well-located  from  openly  available  teleseismic  data. 
These  events  presumably  include  a  few  earthquakes,  and  a  few  tens  of 
explosions  (chemical  or  nuclear).  We  have  no  capability  with  the  data  we 
have  used,  to  distinguish  between  chemical  and  nuclear  explosions.  A 
useful  guide  to  identification  is  the  number  of  each  type  of  event  reported 
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by  various  sources  (see  Table  10),  but,  though  of  some  interest,  it  would 
not  be  prudent  to  rely  too  greatly  on  these  numbers;  nor  can  such  a  guide 
be  used  in  future  CTBT  monitoring.  Our  event  list  can  be  made  the  basis 
for  numerous  projects  to  evaluate  teleseismic  detection  capability,  and  also 
to  evaluate  regional  discriminants.  We  have  accumulated  and  used  several 
hundred  analog  seismograms  for  these  events.  Many  of  these 
seismograms,  particularly  those  obtained  on  ChISS  instruments,  can  be 
used  for  spectral  analysis;  and  we  believe  digital  data  can  also  be  found, 
now  that  we  have  a  list  of  events  as  the  basis  for  data  requests. 
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28/  279;  470; 
■9/  'l29  F‘259;' 


16/  107;  232; 


6.063 


hrmin  jinformation 


Tzero,  s 
i^F^' 


Tzero,  s 

. . 


lat/long 

. " 


mag 


mag  ;  mag 


81-JUI-05 
81 -Sep-36 

81- N0V-19 

82- j"un-1l" 
82-'jui-12 
82-’Sep-64' 

82- 'Sep-T5 

83- Feb-i3 
83^JuF28" 

84- jun-23 

84- ’Augl2® 

85- Jun-27 


85-Jul-ll 

87-Jun^'9 


87- Sep-1 6 

88- Sep;26 

88- Dec-28 

89- Oct-20‘ 


03:59 

"l'2:55" 

05:57 

"10759’ 

"iO:'29' 

05747' 

04:33 

"d3:’02" 

•■q2:57" 


V.Kh.;  S&R 


iV.Kh.;  S&R 
jyxK-s&R 
iRing. 

"]v;kh;7s&R 
lvXh;;S&R 
Tv^kh^-'s&R 
Iv.kh.;  S&R 

'  Iv!  An . 

■jVXh^rs&R 

'IVXhirs&R 


(50,  78) 
.:(50,  .78) 


..(50, ...78).. 
..(5.P.,....78).. 


..(50,.. ..78).. 
...(5.P.,....7.8).. 


(49,  78) 
.............. 


07:30 

■■07:45' 

05:28 

■■i'3:23 


iV.Kh.;  S&R 
lv.Kh.;Sult:K 
iV.Kfi.;  Ring.(96) 

WM . 


5.P.,....7.8 


mpZzsi; 


Station  code 

Name 

Lat 

CHGU 


ALB 


ART 

BAY 

BRV 

CHK 

bza 


KKB 

WR 


KRD 

ms 

MAK 

SVE 


VED 

VOS 

^Ni" 


Distance  j  Azimuth  i  Permanent 
km  I  “from  N  i  or 

temporary 


9.95; 

^;49[ 


Semipalatinsk  |  50  ip 

Skaiistaya  |  49.57 

Ust-Kamenogorsk  i  49 .831 

Ust-kan  I . M;.6p 

Ust-Eiegest . ] . 51 .34 

.iSWiON  JO 

iNovosibirsk  54.85 


88.35 

■■86.'27 


87.19 


81.08 


81.22; 


82.28! 


Alma-Ata 


BoomskoYe 

BaYan-koi 


43.27 


43.03;  80.46 


K.r3?n9.9.or.K3 

kurmenty 


Markatau 


Namangan  ]  4T.pp 


Ortomerke 

42.95 

Podgornaya 

43.33 

Przhevalsk  i 

42.29 

Rybachye  ; 

42.43 

Tegermen 

43.38 

STATIONS  TO  THE  WEST 


Ala-bota 


Artj 

BaYan^Ayj 

Boroyoye 

Chiraloyo 

Dzhukek 


Kashkarbay 

kar-karaiinsk 


Krasny  Kordon 
Kustanay 
Makinka 
Sverdlovsk 


Vedeneyka 

Vostochnaya 

Zerenda-1 

Zerenda-2 


53.63! 


56.40! 

50.82! 

■‘53"6"6i‘ 


53.09! 

Tg'ssF' 


52.96! 

wFliS;' 

■■52y57;" 

■seFsoj" 


52.63! 

"52F7i^ 

52F93[ 

52.88! 


70.92 


58.60! 

■76F2& 


.72 

.61 


69.07 

YsJe 


69.03 

63.40 

■■yQg2 

‘■‘eo.eo 


69.50; 

70.971 

■'69.'(j'5l 

■'69^1^ 


273 


171 


329 


304! 


488! 


1144 


66! 

T 

)!  84! 

T  ! 

ChISS 


ChlSS 


icJoU. 


/f 


66-Oct-29 


66-NOV-19 


66-Dec-26 


03:58:00 


03:58:00 

17:39:38 


03:57:00 


81-JUI-05 


81- Se]j-30 

i.i:.Noy- 1.9 , 

82- Jun-1  i 


03:59:18 


12:55:10 

65:57:14 

16:59:67 


04:33:19 


location; 
. liong) 


eg . L 

9.85i 


49.64; 


71 -Jan-29 

05:03:00 

71 -Apr-09 

02:33:00 

73-Mar-23 

06:30:00 

73-bec-31 

04:03:00 

75-bct-05 

04:27:00 

76-Mar-20 

04:03:39 

79-May-24 

04:07:00 

79-Sep-14 

07:33:00 

79-Sep-15 

04:07:00 

80-Jul-13 

08:10:00 

50.00: 

56!26i 

56.661 

9.82 

6.T6 


deg  I 
77.94i 


78.05 


49.791  78.02: 


49.79:  78.25 


79.84: 


49.85:  78.79 


78.801  3.5i  6. 

78.7l[  10. 6i  10. 
78^6^”  4.1:  12.9 


r\ 

5 

1 

3 

^  - 

1 

8:  li 

3 


8 


talie.  i 


est.  location:  est.  sigma: 

Number  of  stations  used: 

Date 

Time 

Oat,  llongj  l(lat,  Hong) 

!  #  !  E  I  S  !  W  !  N 

deg  i  deg  I  km  I  km 

64-Jun-06 

00:00:00 

49.801(77.40)1  7.81  100. C 

6!  -!  6!  -! 

64-Auq-18  I 

08:00:66 

49.81i(79.10)i  8.5i  100. C 

l!  2i  -!  2i  -! 

65-Feb-04  j 

06:00:00 

49.801(78.90)!  6.4!  40.C 

1!  4!  -!  4!  -! 

87-Jun-29  I 

04:55:08 

50.10!  78.96  !  18.0!  16.C 

li . 2[ . 2l  -!  -! 

Date 

Time 

Station 

AD  ! 

direction 

61-Jun-05 

03:50:00  ! 

MIX 

avr-4  i 

3  i 

North 

69-Apr-64 

03:57:00  ! 

7LG 

avr-3  ! 

6  ! 

North 

74-Sep-27 

07:34:00  ! 

TLG 

avr-3 

30 

North 

81 -May-28 

04:08: 

UKN 

avr-2 

13 

West 

85-Jun-27 

11:57:00  ! 

UKN 

avr-2 

10  i 

West 

89-Oct-20 

13:22:45  ! 

BRV 

avr-2 

21 

West 

iest.  location: 

est.  sigma: 

Mumber  of  stations  used:  1 

Date 

time  !(lat,  !lonq) 

(lat,  llong) 

#  ! 

E  ! 

S 

W  !  N 

deg  deg 

km  !  km  ! 

66-Jun-03 

09:21:10  !  51.68  175.10 

6.5!  13.0! 

3i 

1! 

2 

67-Jul-16 

4:07  !  51.80!  78.90 

25.0!  3.0! 

1 

81 -Mar-31 

07:51:30  !  47.97  !(80. 00) 

10.0!  50.0! 

5! 

-j 

5 

84-Aug-26 

03:32:57  !  50.38  171.64 

2! 

. 21 

- 

-I 

same  event 

alt.  location: !  55.60  !  95.70 

-j  -i 

2! 

21 

-i 

80-Oct-26 

11:49:44  far  from  STS 

(AD  =  813  km)  i 

1! 

.! 

1 

-1 

83-Feb-13 

03:02:09  also  far:  (AD  =  277  or  340  km) 

1! 

. 

- 

-i 

ZO-May:?? 
71 -Jan-29 
71-Apr:09”; 


80-Sep:2p 

sp^Sep^SO" 

80-Sep-36 

8O-N0V-O6' 

8i‘-Mar-31 


81:May-:28 

Si-Jun-OS 

81-Jui-b5 


82-Jun-11 

82-Jui-i2 


82-Sep-04 

82- Sep-i5 

83- Feb-13* 


Date 

Time 

61-Jun-05 

03:50:00  1 

61 -Oct- 11 

07:40:00  i 

62-Feb-02 

08:00:00  I 

64-Jun-06 

00:00:00  i 

64-Aug-18 

06:0( 

65-Feb-04 

65-Mar-27 

06:30:00  j 

65-Jul-29 

03:05:00  I 

65-Oct-14 

04:00:00  i 

66-Oct-29 

03:58:00  i 

66-N0V-I9 

03:58:00  i 

66-Dec-26 

17:39:38  i 

67-Jun-03 

09:21 

:10  I 

67-Sep-02 

04:04:00  i 

68-Oct-21 

03:52:00  j 

68-Oct-29 

03:54:00  i 

68-N0V-I2 

07:30:00  i 

69-Apr-04 

03:57:00  I 

69-Apr-13 

04:04:00  j 

69-NOV-27 

05:02:00  | 

04:03:00 

05:03:00 

02:33:00 


73-Mar-23 

:  06:30:00 

73-Dec-31 

i  04:03:00 

74-Sep-27 

j  07:34:00 

75-Oct-05 

i  04:27:00 

76-Aug-04 

1  02:57:00 

77-N0V-27 

1  03:57:00 

78-Jul-31 

i  08:00:00 

10:40:01 
65;  57:1*2 
05:57:^7 
'l'7:*42:*58 
07:5T:'30 


04:08:28 

03:22ri'8" 

63:59:14 


10:59:07 

16:29:1*8 


05:47:17 
■64*:33:i9 
■‘6*3:62:69 


11.31  4.44 


NORSARlhPS  !ISC 


8.5i  3.21 


81-JUI-05 

i03:59;14 

49.87i  78.99i 

15! 

18! 

3.88!Bal 

0.90!UNE 

0.50!orChE 

81 -Sep-30 

112:55:10 

49.94:  78.901 

1  0! 

15! 

3.96!Bal 

1 .00!ChE 

0.90!orUNE 

81 -Nov-1 9 

i05:57:15 

50.1 1i  78.95! 

12! 

15! 

3.69!N  from  Bal  !ChE 

0.70:orUNE 

82-Jun-11 

h0:59:06 

49.93!  78.50! 

12: 

20! 

4.00!between  B&DiChE 

0.90!orUNE 

82-Jul-12 

i10:29:18 

49.90!  77.90! 

20! 

40! 

4.1 5!N  from  Deg  !ChE 

0.90!orUNE 

82-Sep-64 

105:47:17 

50.06!  78.56: 

10| 

10! 

3.64!between  B&DiChE 

0.90!orUNE 

82-Sep-15 

|04:33:17 

49.851  78.85! 

20! 

40! 

4.18!Bal 

0.50!ChE 

0.80!orUNE 

83-Jul-28 

i03:41:25 

50.07!  78.60! 

12! 

12! 

4.1 7!between  B&b!ChE 

0.90!orUNE 

84-Jun-23 

102:57:16 

49.92!  78.93! 

6! 

10! 

4.34!Bal 

I.OOiUNE 

0.50|orChE 

85-Jun-27 

111:57:04 

one  !station! 

-! 

-! 

3.21!Deg 

0.20!  UNE 

0.20!orChE 

85-Jul-11 

:02:57:02 

49.78!  77.90! 

To! 

3. 94!  Deg 

0.50!UNE 

0.80!orChE 

87-Jun-29 

104:55:11.5 

one  istation! 

-i 

3.21!Deg 

0.40!UNE 

0.70:orChE 

87-Sep-16 

107:30:01 

49.86!  78.73! 

. l'^ 

4.06!Bal 

0.50!  UNE 

0.70!orChE 

88-Sep-26 

i07:45:01 

50.08:  78.80! 

12: 

15! 

3.91!N  from  Bal  iChE 

0.50!EQ  6.4 

88-Dec-28 

105:28:09.5 

49.80I  78.06! 

. 4] 

. 5i 

4.16!Deg 

I.OOIUNE 

!  1 .00! 

89-Oct-20 

113:22:45 

one  !station! 

3.63!Bal 

0.30!Collai: 

)se? 

1 1 


Number  per  year,  of  underground  nuclear  explosions  at  Semipalatinsk 


77.0  77.5  78.0  og  78.5  79.0  79.5 


The  boundaries  of  the  Semipalatinsk  Test  Site,  East  Kazakhstan 


U  A 


ChISS  amplitude  response 


2i}-32U4p-6.4 


ChISS  impulse  response 


4.0  -  6.4  Hz 


3.2  Hz 


1  second 


1.0  -  K6  Hz 


0.8  Hz 


Iih2 


0,6S 


1966-Dec-26  (49.40°,  78.75°) 

(8  km,  15  km) 


1967-Jun-03  (51.68°,  75.10°) 

(6.5  km,  13  km) 


1980-NOV-06  (50.15°,  78.80°) 

(2  km,  4  km) 


l98  1-Jun-05  (49.80°,  /8.»0''j 

(7.9  km,  15.5  km) 


1981-Mar-31 


(47.97°,  80.00°) 
(10  km,  50  km) 


1981-Jul-05 


(49.85°,  78.90°) 
(6.3  km,  10.5  km) 


1981-NOV-19  (50.20°,  78.71°) 

(10  km,  10  km) 


1984-Jun-23 


(49.93°,  79.01°) 
(2.2  km,  4.5  km) 


1987-Sep-16  (49.85°, 

78.79°) 

(6.7  km, 

12.7  km) 

1988-Sep-26 

(50.10°,  78.96°) 

(18  km,  16  km) 

mb(AWE) 


Comparison  of  mb(AWE)  with  mb(ISC)  for  100  Balapan  UNEs 


mb(ISC) 


Comparison  of  mb(Lg)  at  NORSAR  with  mb(P)  from  AWE 


8 


9 


10 


(AWRE) 


m(AWRE)  =  -1  .f  7m(CHISS)  -  i  .02 


282  known  &  located  presumed  UNEs  at  the  Semipalatinsk  Test  Site 


known  &  located  earthquake  (76-Mar-20,  cross);  1  chemical  explosion  (84-Sep-15,  square) 


43  previously  unknown  or  unlocated  seismic  events  on/near  the  Semipalatinsk  Test  Site 


18  UNEs  (circles);  22  ChEs  (squares);  2  UNEs  or  ChEs  (triangles);  and  1  Eq  (cross) 


